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A simple model based on QED is presented for the estimation of contribution of the excited level few-photon 
decays to the radiation escape from the matter in the epoch of the cosmological hydrogen recombination. It is 
shown that apart from the widely studied two-photon decays, some specific 3-photon decays can contribute on 
the level of 0.1% accuracy, required by the recent astrophysical observations. 



I. INTRODUCTION 

Theory of the cosmological hydrogen recombination became one of the most intensively discussed fundamental physics 
problems in the last decade. The interest comes from the accurate measurements of the asymmetry in the temperature and 
polarization distribution of the Cosmic Microwave Background (CMB) 1 1 1, |2|. The launching of the Planck Surveuor which 
enables to perform these measurements with accuracy 0.1% makes the situation even more intriguing. 

The modern theory of the cosmological hydrogen recombination starts from the papers by Zel'dovich, Kurt and Sunyaev 121 
and by Peebles |4|. It was argued that the bound-bound one-photon transitions from the upper levels to the lower ones did not 
permit the hydrogen atoms to reach their ground states: each photon released in such transition in one atom was immediately 
absorbed by another atom. These reabsorption processes did not allow the radiation to "escape" the interaction with the matter 
As it was first established in |[3l, H the two-photon 2s- Is transition presents the main channel for the radiation "escape" and 
formation of the CMB. This transition also led to the final hydrogen recombination. Hence, the recent properties of the CMB 
are essentially defined by the two-photon processes during the cosmological recombination epoch. 

In Is) the importance of the two-photon decays from excited states with n > 2 for the detailed analysis of the properties of 
CMB was noted. Over the past few years the theory of cosmological recombination was essentially detalized by many authors. In 
particular, in |5|, |6| it was demonstrated that the two-photon transitions ns ls{n > 2) and nd — Is can also give a sizeable 
contribution to the radiation "escape". There is a difference between the decay of ns{n > 2), nd states and the decay of the 2s 
state. This difference is due to the presence of cascade transitions as the dominant decay channels in the cases of ns{n > 2) and 
nd levels. For the 2s level the cascade transitions are absent. The cascade photons can be effectively reabsorbed and therefore 
the problem of separation of the "pure" two-photon emission from the cascade photons arises in connection with the "escape" 
probability. This problem was intensively discussed during the last decade |7 1-| 12 1. As it was proved in 1 1 1 1 the separation of the 
"pure" two-photon emission for the ns{n > 2) and nd levels is an ambiguous procedure. First this ambiguity was established for 
the two-photon transitions with cascades in the highly charged ions 1 13|. To reach the level of accuracy 0.1% for the theoretical 
description of the properties of the CMB many effects should be taken into account in the astrophysical equations describing the 
radiation "escape" process: consequences of the universe expansion, thermodynamical properties, induced radiation, processes 
of the electron scattering, Raman scattering etc. Detailed analysis of the various distortions of the resonant optical line spectra 
is also required 1 14|, including the nonresonant corrections 1 15 1. 

This very complicated construction requires a careful treatment of the basic principles which this construction is standing 
upon; these principles are given by the Quantum Electrodynamics (QED). In this paper we will analyse these principles and 
demonstrate that following them one can find some additional effects, small but probably sizeable at the level 0.1%. Our 
treatment will remain in the frames of QED applied to the free atoms in the field of photons; the astrophysical aspects will 
be restricted to the introduction of temperature (i.e. thermodynamical equilibrium). Actually we consider the model universe 
containing two atoms only. The first atom is in excited state and the second one in the ground state. The first atom emits the 
radiation and aiTives in the ground state too. If this radiation is not absorbed by the second atom, both atoms appear to be in 
the ground state: recombination occurs and the radiation has "escaped" the interaction with the matter This "escape" does not 
coincide with the definition adopted in astrophysics and has sense only within the frames of our model. We assume however that 
our model can reproduce correctly the relative role of the higher excited states compared to 2s state in the process of hydrogen 
recombination. 

Our paper is organized as follows. In Section II we describe the process of the photon scattering (i.e. the process of scattering 
of the photons emitted by one atom, on another atom) in QED and apply this description to the rescattering of the Lyman- 
alpha photons. In Section III we investigate how the two-photon emission from the ns, nd levels is absorbed in the one-photon 
transitions. This investigation will enable us to compare the probability of the radiation "escape" from the ns{n > 2), nd levels 
with the "escape" from 2s level. Section IV is devoted to the studies of the multi-photon (i.e. 3-, 4-photon) transitions in the 
two-photon approximation. We will show that these transitions can give a non-negligible contribution to the radiation "escape". 
Section V contains discussion of the results and conclusions. In Appendix A we give a rigorous QED derivation of the Lorentz 
contour for the one-photon transition between the two arbitrary excited levels; such a derivation was yet absent in the literature. 
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Appendix B is devoted to the derivation of the basic formula employed in section III. 



II. QED THEORY FOR THE PHOTON RESCATTERING ON AN ATOM 

A. Emission line profile for the transitions between two arbitrary levels 

A quantum mechanical phenomenological description of the line profile (Lorentz profile) is known since 1930 lfT6l . ifTTl . A 
QED derivation of the Lorentz profile for the transition between two excited states was given in ifTSl . ifTSl . The corresponding 
expression for the transition a ai looks like ifTSl : 

^T/f/ f, ,\ — ^ ^aai^aiao(ra + ^ai) 

27r iaiai (i^ - UJaaiV + li^a +^air 

In Eq. ([T|i it is assumed that the lower level ai decays in turn to the ground state aq via one-photon decay. Here Fa, Ta-^ are 
the total widths of the levels a, ai and Taai, ^aiao the partial widths corresponding to the transitions a — J' Oi and Oi — )■ Gq. 
Emission probabilities are connected with the partial widths via equalities 

aa\ — Tqq^ , ^^aiao — Ta^aQ . (2) 

Finally, ujaa^ = Ea + La — — Lai, where Ea, Ea^ are the one-electron energies and La, La^ represent the Lamb shifts of 
the levels a and ai. Thus, in principle, the line profile for the transition a — !• ai depends on the further decay channel for the 
lower state ai. Actually this is the dependence on the branching ratios baai ~ ^aai/^a and &aian = ^aiao/^ai- In the simplest 
case when both levels a and oi have only one decay channel baai — baiao — 1 EQ- fl} simplifies to 

/ ^ 1 {Ta+Tajduj 

2tT (uJ - UJaai V + i[^a+ ^ai Y 

Here the dependence on the state ao disappeared totally. In Appendix A we present also the derivation of the Lorentz profile for 
the most general case when the lower level ai decays not directly to the ground state ao but to the intermediate state 02, then 
to the lower intermediate state 03 and so on. The total chain of decays (the cascade) is a ai ai — J- 03 ^ ... qq. It is 
assumed that these decays are of one-photon type. 

An important question is: how far from the resonance the wings of the Lorentz profile can be extended? The answer depends 
on importance of the so called nonresonant (NR) corrections which distort, in principle, the line profile. The NR corrections were 
first introduced in 1 19| and recently discussed in connection with atomic laboratory experiments in 1201 - 1221 . In the astrophysical 
aspect the role of the NR corrections was studied in U5J . According to these studies, the Lorentz profile can be extended far 
from the resonance (actually, to infinity) without any serious errors. We will use this extension throughout this paper. 

Employing the extension of the profiles discussed above we choose the normalization condition for the Lorentz profile as: 

dW{uj) = L{Lu)duj , (4) 



^ LaaA^)'^ = 1 ' (5) 


in case of Eq. (|3]l and 

Laai{ao)i^)di^ ^ baaiba^ao , (6) 







in case of Eq. ([T]|. Eqs (jSj), (|6| represent the absolute probability for the photon to be emitted via the transition a — ?> ai with any 
frequency value. If there are no other decay channels apart from a — > ai, this probability equals to L If such decay channels 
exist both for a and ai levels, this probability is defined by the product of the branching ratios fcaai^aiao- 



B. Reemission of the photons emitted by one atom by another atom in the same transition. 



Now we assume that the radiation emitted in the transition a ^ ai and having the frequency distribution defined by Eq. ([3]l, 
is absorbed by another atom via transition ai a. The simplified form of the Lorentz profile takes place, in particular, for the 
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most important for the cosmological recombination Lyman-alpha line. The absorption Lorentz profile is defined by the Eq. ([3]) 
and the absorption probabilities are connected with partial widths via 



9ai 
9a 



where ga^, ga are the degeneracies for the states ai, a. If the absorbed photons originate from the emission line a 
another atom, the probability of the absorption and hence reemission of these photons should be defined as 



(7) 
ai of 



aai 



I aai 



(8) 



where laai {^) = i^a + Tai )Laai {^) IS the dimensionless distribution of the incident photons. This function is normalized 
according to the condition 



oo 

J laai iuj)du} = Ta + Ta^ 



(9) 



The formal proof of Eq. ([8]) on the basis of QED is given in Appendix B. 

The frequency distribution of the emitted photons was first introduced in the QED S-matrix theory in |fT9ll and later employed in 
1, ll24l for the studies of the multiple photon scattering on the hydrogen atom. 
In this way we can define also the probability of the photon emission after the multiple (n-fold) scattering: 



(10) 



For n = 1 the expression ( 10 1 reduces to Eq. (j5]l. The integral can be extended to the interval — oo < cj < +oo since the main 
contribution comes from the pole in the expressions ([T]i or ([3]l for Laai ■ Then we can evaluate the integral in the complex plane. 
Employing the equality 

1 1" 1 

(11) 



and using the Cauchy's formula we find 



[iU - Cjaa, - 5 (Fa + )] " [C^ - cZ-aai + | (Fa + T,, )] ' 



{n-l)\ 



(n-1) 



.i + i(r„+r„,) 



2tt 



where f...]*^" denotes the {n — l)-fold derivative with respect to the variable to. The evaluation in Eq. ( 12 1 results 



Wn)^ (2n-2)! 1 
""^ ((n- 1)!)2 (27r)"-i' 



(12) 



(13) 



For all n > 1 it is easy to check that Xqq^ < 1, so that we can interprete x'^^ as the absolute probability of the absorption in 
the process of rescattering and reemission. Then the quantity 



((n- 1)!)2 (27r)"-i 
can be interpreted as the probability of the radiation "escape". In particular, 

1 



y(2) = 1 



0.682 



TT 
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(14) 



(15) 



(16) 



So, in our simple model the probability of the radiation "escape" directly via the arbitrary one-photon transition is 0.682 already 
after the first rescattering and becomes close to 1 with the increase of number of rescatterings. This result does not depend on 
the particular transition and corresponds also to the Lyman-alpha transition. We should stress that these estimates cannot replace 
the accurate astrophysical approach to the problem of the photon rescattering on the matter and are presented here only to make 
the further derivations more obvious. 



4 



III. RADIATION "ESCAPE" IN THE TWO-PHOTON TRANSITIONS 



A. QED theory of the two-photon transitions 



Quantum mechanical theory for the two-photon transitions was first developed by Goppert-Mayer f25l and the first evaluation 
of the two-photon 2s — > ls+2'^{El) decay rate in hydrogen was performed by Breit and Teller |26 1. The accurate nonrelativistic 
calculation for this transition rate was given in ll27l . The fully relativistic calculations valid also for the H-like ions with arbitrary 
nuclear charge 1 < Z < 100 were performed in Il28l - ll30l . The most accurate recent calculation with QED corrections one can 
find in [i311. 

The modifications of the theory necessary to describe the two-photon transitions with cascades were discussed in 1131 (see 
also LiU, Li8J). 

The transition rate for the 2s Is + 2j{El) transition in H atom looks like (in atomic units) 



ZlTT 



(17) 



where 



o / \ _ < Ris\r\Rn'p \r\R2s > 

n p ^ 



(18) 



<Rn'i'\r\Rni >= / r^Rn'i'{r)Rni{r)dr 



(19) 



cjq = £^2s ^ Els, Rni{f) is the radial part of the nonrelativistic hydrogen wave function, Eni are the electron energies for the 
hydrogen atom and a is the fine structure constant. Due to the absence of the energy levels between 2s and Is (i.e. cascades) 



the denominator in Eq. ( 18 i has no zeros. 



The total decay rate for the two-photon transition 2s — > Is can be obtained by integration of Eq. ( 17 1 over the entire frequency 
interval 



W2s,is = \ I dW2s^s = 8.229 s-i 



(20) 



An expression for the transition rate W^s.is in the presence of cascades was given in fTTl : 



(21) 



where 



jjrT{cascade) 



4 r3, + r2p 



277r r 



2p 



(II) 



{Rzs{r)\r\R2p{r)){R2^{r')\r'\RUr')) 



E2p — Ey,s + CJ — I (F; 



2 l.^ 3s ^ 2p) 



duj - 



+ 



27n 



(IV) 



{R3s{r)\r\R2p{r)){R2p{r')\r'\RUr')) 



E2p — Els — UJ — |r2p 



(22) 



(II) 

277r J 

(IV) 



did 



duj 



+ 



27t: 



iO^iujQ - Ujf |S'is;3s(w) + Sis:3sii^a - duj 



(I+III+V) 



(23) 
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_l„r{interference) _ f 4:01^(010 — Lpy 
""3s;ls 



/ 



2777 



Re 



(11) 

2Tjt 



{R3s{r)\r\R2p{2r)){R2p{r')\r'\Ru{r')) 

E2p — E^s + LO — |r2p 



Re 



(IV) 



{R3sir)\r\R2A2r)){R2Ar')\r'\Ris{r')) 
E2p — Els — OJ — |r2p 



Sis-Li^) + Sis-M^o - uj) 



duj 



duj (24) 



Here S'{^^j^(a;) is the expression (18i with the term in the n'p — 2p excluded from the summation, luq = E^g — Eig. The 
intervals of integration over w (I)-( V ) are defined as 



(I) o<w<woi-^(r2p + r3,) , 



(25) 



(II) Woi - l{^2p + Tss) < W < Wol + ^(r2p + Tss) 



(III) Wol + l{T2p + T^s) < W < Wo2 - ^(r2p + ^3s) , 



(IV) UJ02 ~ l{T2p + Tss) < W < W02 + l{T2p + Tss) 



(26) 
(27) 
(28) 



(V) a;o2 + ^(r2p + Tss) < w < Wo 



(29) 



where wqi = E^s — E2p, u!q2 = £'2p — £^is are the frequencies for the two links of the cascade, I is integer chosen to separate 
the cascade contribution from the "pure" two-photon contribution. As it was shown in ifTTI . this separation is not unique, i.e. 



the contributions W^^'^i^"''^'^' 
^Ss'is remains invariant: 



, W, 



{pure2'y] 



3s, Is 



, w. 



'inter f erence) 



3s, Is 



change essentially depending on the choice of the I, but the total sum 



(27) 



3s,ls 



dW3s isduj ^ 6.317 • 10'' 



(30) 



Note that the factor '"^p^'"^'' in the first line in Eq. (22) was omitted in ifTTl . so the numerical value Eq. (30) also was different. 
We should stress, that in our derivations the total width of the level T^s does not coincide with the value given by the Eq. (30), but 
coincides with the transition rate v'^^'^lp^as in the atomic spectroscopy. The value r3^^2p' principle, defines the total width in 
the laboratory experiments when the one-photon transition rate for the photons with frequency uJ3s,2p = E^s — £'2p is measured. 
The one-photon decay 3s — ^ 2p appears to be faster than the decay 3s ^ 2p ^ Is due to the destructive interference of the 
cascade decay with the "pure" two-photon decay in Eq. (30). However this difference can be traced only in the fifth digit. The 
same picture holds for the decays of the other ns{n > 2), nd levels. Similar expressions can be written for the transition 3c?— Is 
with the cascade 3d — 2p — Is, for the transition 4s — Is with two cascades 4s — 3p — Is and As — 2p — Is, etc. 



B. Radiation "escape" via two plioton decays 
We define the "escape" probability for the incoming two-photon radiation via the Lyman-alpha channel similarly to Eq. (j8]l: 

^Sls l2s,ML2p,ls{uj)duJ , (31) 



where /2s,is(w) = dW2s,is{oj), ujq — i?2s ~ Eis- The result of the integration 

X^fJ] = 6.50 • 10^22 (32) 

shows that the two-photon 2s — Is radiation emitted by one atom cannot be absorbed by another atom. This means that the 
radiation "escape" via the two-photon 2s — Is transition is absolutely full: 

^(2)27 _ 1^(2)27 _ 

^2s,ls - ^2s,ls - ^ • (-^J) 
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The superscript (2) here, as in Section 11, means that we consider only one scattering (reemission) of the photons. This is enough 
to understand the relative importance of different decay channels for the radiation "escape". 

Now we can repeat the same for the transition 3.s — > Is + 2'^{El). In this case we evaluate the probability of reemission of 
the 3s — Is two-photon radiation via all possible one-photon transitions within the frequency range [0, ojq], i.e. 3s — 2p, 3d — 2s, 
3p — 2s and 2p — Is 

^H]u = ^ J hs,ls{uj) [^3s,2p(w) + i2p,ls(w) + i3d,2p(w) + i3p,2s(w)] duJ , (34) 


where ujq = E^g — E'l^.The numerical result is 

X^^Jis = 0.00497 . (35) 

The value X^^^'^] is much larger than xS^^J^] but still essentially smaller than 1. This means that the "escape" probability is very 
high; 

Y^^^^^ = 0.99504 . (36) 

The same picture holds for the two-photon decays of the other ns {n > 2), nd levels: for the transition 3(i — > Is + 2^{E\) 
which occurs as one cascade with two links; for the transition 4s Is + 2'y{El) which includes two cascades each with two 
links As ^ 3p ^ Is and 4s — > 2j3 — > Is and for the transition Ad ^ Is + 2^{E1) which includes two cascades each with two 

links 4(i — > 3p ^ Is and Ad ^ 2p ^ Is. The corresponding total decay rates W3^^]^'s, wi^^^,, W^4^'^is and the total widths of 

the levels T^d^ r4s, r4c; as well as the probabilities of the reemission ^^^^f^ "escape" probabihties Y^'^^^ are given in 

Table 1. 

Thus, all the levels ns, nd with n = 3, 4 seem to be nearly as effective for the radiation "escape" via the two-photon transitions 
to the Is state, as the 2s level. The smallest "escape" (difference about 3% with 2s level) occurs for the 3d level. This corresponds 
to the maximum "death probability" for the Lyman-alpha photons due to the transitions to 3d level, as found in L14J . However, 
the role of all these levels in the cosmological radiation "escape" is strongly suppressed by the thermodynamical factor (see 
Section V). 

Two comments are necessary concerning the accuracy of the results given above. First, we have fuUy neglected the two- 
photon transitions other than ElEl. For the neutral hydrogen atom it is well justified (see, for example, ll32l ). Second, we have 
neglected the difference between Vaa' and in Eq. ([T]) for the one-photon transitions 3s — 2p, 3d — 2p, As — 3p, As — 2p, 
Ad — 3p, Ad — 2p. For example, in Eq. ( 34 1 the Lorentz profile L-^sap should be defined as 



-, r(i7) p(i7) I r(i7) 



(17) 



As we have discussed above in our derivations we have to put Fss = ^zs ip- In the latter equality the transitions 3s 
Is + j{Ml), 3s 2s + 7(Afl) are neglected, which give only extremely small contribution to the total width T^s 1331 . 



IV. RADIATION "ESCAPE" IN THE MULTIPHOTON TRANSITIONS 



A. Contribution of the 3-photon transitions 

In ll34l . II35I it was suggested that the multiphoton transitions which contain cascades with the two-photon links can also 
contribute to the radiation "escape" in the process of the cosmological recombination. This approach was called "two-photon 
approximation" since the contribution of the "pure" multiphoton transitions with the number of photons more than two were 
neglected. One of the examples described in 1341 was the two-photon approximation for the 3-photon 3p — > Is transition. The 
3p level decay can occur either as a one-photon transition 3p — > Is + 'y{El) or as a 3-photon transition 3p — > Is + 37(i?l). 
These channels do not interfere due to the different number of photons in the final state. The one-photon decay rate is 

W^l'^l = l95Mma^{aZ)^ r.u. = 1.67342 • 10*^ s~^ , (38) 

where m is the electron mass, a is the fine structure constant, Z is the charge of the nucleus (Z = 1 for the hydrogen). 
The 3-photon decay rate 3p — > Is + 37(i?l) consists of the "pure" 3-photon contribution, two cascade contributions 3p 
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2s + 7(i?l) ^ Is + 2^{E1), 3p 2p + 2^{E1) -)■ Is + j{El) and the interference terms. The "pure" 3-photon contribution 
to the decay rate is of the order ma^{aZ)^ r.u. Three-photon contribution was evaluated in |36| for the 2p Is + 3'f{El) 
transition which does not contain any cascade contributions: 

Mlss = 0AM6ma^{aZ)^ r.u. (39) 

In principle, for 3p Is + 3j{El) transition rate the contributions of the "pure" 3-photon decay channel and the cascade 
contributions are again inseparable, similar to the case of the 2-photon 3s — Is + 2j{El), — Is + 27(i?l) transitions as 
discussed in Section III. However, unlike the two-photon decays in Section III, where at the level of accuracy of the "two-photon 
approximation" we were interested in the all contributions, in case of the 3-photon transitions at the same level of accuracy we 
have to keep only the cascade contributions and neglect fully the "pure" 3-photon contributions and the interference terms. This 
simplifies our task and retaining only the cascade terms, we find [34J 

w(37) _ yr.(27) , ^^3p,2s (27) ..qs 
^ 3p 

Note, that in ll36l the right-hand side of the equation corresponding to Eq. (40) contained the wrong factor |, the same concerns 
the equation for the other two-photon decay rates. This mistake was noticed in (37]. The order of the magnitude of the 3-photon 
cascade transition rate is defined by the fact that the cascade transition rate is determined by transition rate of the slowest cascade 
link, i.e. in our case by the two-photon transitions. 

The total decay rate of the 3p level T^p is defined as (see discussion concerning the width T^s in Section II) 

^3p = 4l% + rg%. (41) 

The two-photon transition rate VFgp^p should be evaluated similarly to Wj^^^^^^ transition rate since it is "pure" two-photon 
transition rate. Hence 

Wtp!i = 1 1 dW^'pHi^) , (42) 



where ujq = E^p — E2p. The two-photon frequency distribution dVFjp^p = /3p.2p looks like (in a.u.) 



97r * 152 



^ duj , (43) 



where 



^3pM^)-[^^lg'0 0^1<lh,m,,J2^ E,^, ^ - E^p + CO ^ 

n' s ^ 

_L/'o /^■^/-''^""'■■'p ri2m,„ ^ < R:ip\r\Rn-d >< Rn'd\r\R2p > ... 



^ En-i- Ej,p^uj 

mi^ n'd 



00 

<Rn-v\r\Rni >= j r^Rn'i'{r)Rni{r)dr , 



(45) 



I mi 



h mil '2 ni(; 



are the Clebsh-Gordan coefficients, mi are the angular momentum projections. 



Unlike Eq.(17i there are no resonant terms with zero denominators in Eq.(44l. Therefore we can perform the summa 



tion over n's, n'd in Eq.(44 1 explicitly using the known expressions for the Coulomb Green function l,38J . Performing also the 
summation over all angular momentum projections in Eq. (|44]l, we arrive at 



dW. 



(27) _ 8a;^(a;o - a;)' 



3p,2p 



^5^^^ {2^ll{v) + 76/|(z.) + 2Ml{u') + imh{u')l2{u') + 76/|(i.') + 

+mio{u){m2{u) + + QhW)) + 12/2 (i^) (5/0 (i'') + 21/2(i^'))) , 



(46) 
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The integrals Iq, I2 in Eq.(46 1 are defined as 

loH - 
I2M 





00 00 



dridr2r^r'^R2i{ri)go{v]ri,r2)Rzi{r2) , 
dridr2r^r^R2iiri)g2{v;ri,r2)R3i{r2) , 



(47) 







Z/ -y/ i?3p — u!, v' = Z/ \J Ej,y — bj' , uj' = ujo — oj and the radial part gi{i^; r, r') of the Coulomb Green function is 



AZ f 'i 



-T^rr exp 



V I ^ I2l + l + n)\in + l + l-u) 

^ n— 



(48) 



In Eq. (48 1 i^/^^ are the Laguerre polynomials The radiation "escape" via the 3-photon transition 3p — ^ Is + 37(i?l) should 
be defined as 



(49) 



^(2)37 _ \ 











'^'^3p,2s 



y(2) 
^2s,ls 



(50) 



In Eq. ( 50 1 the decay rate W, 



(17) 

3p,2s 



r 



(17) 

3p,2s 



2.24603 • 10^ s-\ the decay rate = 1.67342 • 10* s'^ and according 

to Eq. (41 1 r3p = 1.89803 • 10** s^^. Note that the difference between Wg^^^^ and becomes much more significant than 



the difference between w'^^i^, ^isll' ^idll ^zd^ '^^s^ respectively. The direct one-photon decays of the levels "id, 
4s, \d are very much faster than the corresponding cascades (see Table 1). This happens because these cascades contain a "pure 
two-photon" link, which transition rate is as slow as the 2s — Is transition rate. Another situation occurs for the the 3-photon 
decays 4p ^ Is + 37(i?l), 4/ Is + 37(i?l). The formulas, similar to Eq.(40 1, for these transition rates look like 



(37) 



4p,ls 



w, 



(27) 



w. 



(17) 



4p,2p 



3p,ls w(27) 
3p 



^^ipM „^(27) 
^^3,d,ls 

-L 4p 



i 4p 



l^V^ (27) ^4p,2. (27) 

-T^ * * 3s,ls ' ^* 2s, Is 

i 4p i 4p 



(51) 



w(37) _w(27) , W(27) I ^''3p,ls (27) .^2^ 

^ 3p 

The two-photon Hnks of the cascades Ap 2p + 2-i{El) Is + 37(£'l), 4/ 2p + 27(£'l) ^ Is + 37(£:i) are 
similar to 3s Is + 27(£'l) transition, i.e. they represent themselves the cascade transitions: Ap ^ '5d ^ 2p and 4/ ^• 
3d 2p transitions, respectively. Accordingly, these links have to be evaluated in the same way as 3s Is + 27(i?l) 
transition and are comparable by the magnitude. With it again the total decay rates W^i^^p, W^'^^p differ from the total widths 
= rip^i^, + r[l% + r[l% + rf^^]^ and r^f = r[)^^l^, respectively. The two other decay channels of Ap, Af levels, i.e. 
Ap ^ 3p + 2j{El) — > Is + 37(i?l) and 4/ — ^ 3p + 27(£'l) — >■ Is + 37(i?l) contain "pure two-photon" cascade Hnks and, 

therefore, have very low transition rates compared to T^p, T4f. The results for the X^f], vj^^^^ and xj^j'fj, Y^j''^^ are also 
given in Table 1 . 



B. Radiation "escape" in the four-photon transitions 

Now we turn to the 4-photon processes and consider first 4s — > Is + 47(i?l) process. This process can occur parallel to 
the 2-photon process 4s —>^ Is + 27(i?l) process, described in Section III and does not interfere with the latter due to the 
different number of the photons in the final state. As in the case of 3-photon decays we neglect the contribution of the "pure" 
4-photon decays and consider only the cascade contributions with 2-photon hnks, namely 4s 3p + 7(i?l) ^ 2s + 27(i?l) 
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Is + i-yiEl), 4s ^ 3s + 2-f{El) 2p + 3j{El) -)■ Is + 4:j{El) and 4s 3p + 'y{El) 2p + 3'y{El) ~> Is + 4j{El). 
Then the expression similar to Eq. (pOb arises j 



4s i 4s i 3p 



In Eq. (|53| wj^^^Sp = ^4s?3p ^^^^^ °f ^1^^^ '° 



-L 4s — i 4s, 3p + -L 4s, 2p ■ y-'^' 

Similarly, W^l'^l = v'^^j^l and Tgp is defined by Eq. 

The two-photon differential transition rate (frequency distribution) dVl^i^ g^, can be obtained from Eqs |l7| , |l8| with the 
replacement 2s — > 4s, Is — > 3s and the value of the decay rate magnitude can be found in ll35l . Then the radiation escape" via 
the 4-photon transition 4s — t- Is + 4:'f{El) is 

^(2)47 _ _ 1^.(2)47 

^4s,ls — ^ ^4s,ls ' y-'-'' 



1 



^4s!lJ " ^ / ^4s,3s(w) [i4s,3p(^) + i4p,3s(w) + L4/,3d(w) + L4d,3p(^) + i4p,3d(w)] duJ + (56) 







1 w(i'y) 7 w^i-^) w^i-^) 

J- '^'^4s,3p T t \!T ^ I r / \ I r / M J i '^'^4s,3p '^'^3p,2s ^(2)27 
^ / -'3p,2p(w) [£3p,2s(w) + -t-3s,2p(w) +i3d,2p(w)J aw + — — -'^2s,ls ' 



i 4s J i 4s i 3p 



where ujq = E/^s — i?is- Inserting all the necessary frequency distributions for the one-photon transitions in Eq. (56 1, we obtain 
the result, given in Table 1 . 

The results of similar derivations performed for the 4-photon 4c? Is + 47(£'l) transition see also in Table 1. Note that 
for these 4-photon decays as for the 3-photon decay ip — > Is + 37(i?l) the total decay rates are much smaller than the total 
level widths r4s, Vid- With the calculations performed in Sections III, IV all the contributions to the radiation "escape" from 
nl{n < 4) levels in the "two-photon approximation" are exhausted. 



V. THE ROLE OF THE EXCITED LEVELS IN THE RADIATION "ESCAPE" IN THE EPOCH OF THE COSMOLOGICAL 

RECOMBINATION 

In this Section we will estimate the relative importance of the multiphoton (two-, three-, four-photon) decays of the ns, np{2 < 
n < 4), nd, nf{n < 4) levels in the radiation "escape" in the epoch of the cosmological recombination. We assume, as usually 
that the thermodynamical equilibrium existed and the electron temperature was approximately equal to the photon temperature 
IHO). For defining T-y we employ the formula BOl 



T7=To(l 



(57) 



where Tq = 2.725K is the recent CMB temperature and z is the redshift, which for the estimates can be taken as z 1000 for 
the cosmological recombination epoch. 

Our aim is to compare the role of the excited states in the radiation "escape" with the well known role of the 2s level 13], iU. 
For this comparison any characteristics of the radiation "escape" can be employed, which we will not specify here. Denoting 
this characteristics as i?2s we can suggest the following formula for the estimate of the relative role of the excited states 



R — R2s 



1 



(2)27 



(2)27 



5d^(2)27 
„ -'4(i.ls 



3s, Is 



5p^(2)37 
„ ^4p,ls 
i/s 



3d^(2)27 

r, 3rf,ls 

ys 

g,V,ls 



(58) 



In Eq. (58 i the degeneracy numbers for nl states are: — 2, gp = 6, gd = 10, g/ = 14. The "escape" probabilities Yni.is are 
compiled in the Table 1. For the rough estimate we can set all the ^,[^43^ values equal to 1 (see Table 1). For Tg = 2725K the 
first exponential factor in the right-hand side of Eq. ([58|) equals to 



0.00032105 



(59) 
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and the second exponential factor equals to 

g-^rf^ ~ 0.000019225 . (60) 

Hence, 

R = i?2s [1 + 0.0018501 + 0.0002941] = i?2,. [1 + 0.0021442] . (61) 



The numbers in of Eq. (59 1 mean that the radiation "escape" from all the excited levels with n = 3 can contribute at the level 
0.18% and the radiation "escape" from all the levels with n = A can contribute about 0.03%. Both these numbers may be 
essential on the recent level of accuracy of the astrophysical observations. It also evident that the levels with n > 4 can hardly 
give sizeable contribution. From this formula we excluded 3-photon decay of 3p level and 4-photon decays of r4s, levels 
which were considered in subsections A, B of Section IV. The branching ratios for these decays, according to the Table 1, are 
too small. To contribute essentially to the "escape" probability the cascade transition has to contain two-photon link which is 
not "pure" two-photon. The only 3-photon decays which satisfy this condition are the decays of T4p, levels. In other words 
essential contribution comes from the 3-photon cascade processes. The 4-photon decays, compatible with this condition can 
arise only for the nl levels with ti ^ 5, for example 5s Is + 47(£^1). 



VI. CONCLUSION 



In this paper we have presented a simple model, based on QED, to estimate the relative role of the excited level few photon 
decays in the radiation "escape" in the cosmological recombination epoch. These estimates cannot replace the accurate solution 
of the astrophysical balance equations but can give a hint which processes (decays) should be included in the these equations. 
In particular, it appears that the 3-photon cascade decays of the nl{n < 4) levels can give a contribution comparable with the 
widely discussed contribution of the two-photon decays. Our studies are based on the "two-photon approximation" when we 
take into account the cascades which apart from the one-photon links have also one two-photon link. This approximation seems 
to be well justified due to the relative smallness of the "pure" 3-photon, 4-photon etc. processes and due to the smallness of the 
processes with several two-photon cascade links. 

In our investigations we also used an idea of the photon reemission probability first introduced in QED by F. Low |19|. 
According to this idea, if an atom emits a photon which frequency is distributed as a result of the preceding absorption, the total 
probability of the photon emission can be smaller then 1 . The deviation of this probability from the unity reflects the fact that the 
incoming photon was not necessarily absorbed. This deviation we associate with the radiation "escape". We assume also, that 
the use of one re-emission (Y^^^''^^ value) is enough to characterize the importance of a certain decay channel for the radiation 
"escape". 

The total contribution of the excited ns{n > 2), nl{l — 1, 2, 3, 4, n < 4) levels to the radiation "escape" compared to the 
contribution of 2s — Is + 2j{El) process according to our estimates reaches the value 0.21% which is not negligible in view 
of the growing accuracy of the recent astrophysical observations. The smallness of this contribution is due exclusively to the 
relatively low equilibrium temperature during the epoch of the cosmological recombination. 
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TABLE I: Contributions of the multiphoton cascade processes, having one two-photon link, to the radiation "escape". Here k is the number 
of photons; nl - initial state of an atom; W^'^^^ij, - total transition rate for the cascade transition;W^*]j'*, - transition rate via "pure two-photon 
links", r„; - total width of the upper level in the two-photon link; X^^''^^, - reemission probability for the photons from the two-photon link 
of the fc-photon cascade; Y^iJ^li, - "escape" probability for the photons from the two-photon link of the fc-photon cascade. 





k 


nl 


iy<tT',s-i 

rii,ls ' 


nl,ls ' 


r„i,s"^ 


^(2)fe7 
nils 


y(2)fe7 
nils 


1 

2 
3 
4 
5 


2 


2s 
3s 
3d 
4s 
4d 


0.06317 X 10* 
0.64686 X 10** 
0.04171 X 10* 
0.26013 X 10* 


8.22935 


8.22935 
0.06317 X 10* 
0.64686 X 10* 
0.04416 X 10* 
0.27677 X 10* 


6.39353 X 10"^^ 
0.00497 
0.04652 
0.00431 
0.02118 


1.00000 
0.99504 
0.95349 
0.99569 
0.97882 


6 
7 

8 


3 


3p 
4p 
4/ 


0.003929 X 10* 
0.784812 X 10* 


1.01909 


1.89803 X 10* 
0.81311 X 10* 
0.13795 X 10* 


2.33432 X 10"^^ 
0.011005 
0.08435 


1.00000 
0.98995 
0.95348 


9 

10 


4 


4s 
4d 




0.61571 
0.41132 


0.04416 X 10* 
0.27678 X 10* 


2.47954 X 10"^^ 
1.75724 X 10"^^ 


1.00000 
1.00000 
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Appendix A: Resonant scattering of photons on an atomic electron and the Hne profile for the emission process 



We consider first a n-photon elastic scattering process, depicted in the Feynman graph Fig.l. The corresponding S-matrix 
element can be written as 

S{x2„-lX2nh,,2r^A^^2^'''"''Hx2n)'ipao{x2n) , (A.l) 

where S"^^"^ is the 5-matrix of 2n-order, e is the charge of electron, i^aix) = V'a(?06 the solution of the Dirac equation 

for the atomic electron, Ea is the Dirac energy. ■0a (x) = V'Ito is the Dirac conjugated wave function with ^/^J being it's 
Hermitian conjugate, 7^ = (70 > 7) are the Dirac matrices. Wave function of photon Afj^{x) looks like: 

^e, k) ^ ^ e(A)g.(fcf-.t) ^ ^ e|,^)e-'-* A^f' (f) , (A.2) 

where e^j^^ is the photon polarization four-vector, k — {k, uj) - is the photon momentum four-vector (k is the wave vector, 
Lj = \k\ is the photon frequency), x = (r, t) is the coordinates four-vector (r, t are the space and time coordinates). Function 
(A.2 1 corresponds to the absorbed photon and the function A* (x) corresponds to the emitted one. S{xiX2) is the Feynman 
propagator for atomic electron. In the Furry picture, the eigenmode decomposition for this propagator reads: 

00 _ 



2n J uJi — £js(l — lO) 



where summation in (A.3 1 extends over the entire Dirac spectrum of the electron states s in the field of nucleus. Integration over 
frequency and time variables in (A.l 1 leads to 



{Eao - Es„ + ^i)i^ao " £^s„+i + J2iZn+2 ^i) ■ ■ ■ i^a,, - E2n-1 + i^2n) 
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The amplitude of the process of elastic photon scattering is related to the S'-matrix element via 



2n 



\z— 1 z— n+1 

The energy conservation in this process is implemented by the condition 

n 2n 



^ = ^ (A.6) 



i—1 i—n+1 



and the resonance frequencies are given by 



^2n — Wl — Ea^ — Eag 
W2n-1 = ^2 = Ea2 — Ea-^ 



Wn+1 — (^71 — Ea„ — Ea„_-i 

[Ea^ =Ea). (A.7) 



Accordingly, we will obtain for the scattering amplitude the expression 



Tjsc.(2n) _ {Uui)a.QSi{Ui^^)siS2 ■ ■ ■ {Ui^^)sn-iSn 



Si,S2,...S„ 
S„ + l,...,S2ii-l 



(C^Wn+i )s„S„ + i • • • {UuJ2n-l)s2n-2S2n-l {Uu.2js 



(A.8) 



where we abbreviate the transition matrix element as 

{U^)ab={aA^''''^)ab ■ (A.9) 



In the resonance case we have to retain the terms si = ai, S2 — a2, . . . , s„ = a, . . . , S2n-2 = Q,2j S2n-i — o-i in Eq. (A.8 1, 
which yields 

Tjsc.{2n),res _ (^cji )aoai (t^tj2 )"i"2 • ■ ■ (^w„)an-ia 

~ (Eao - Ea, + L^lKEao " + ^1 + C^a) • • • (^ao ^ + ^^^^ iO,) 

{Uu^)aa„-i ■ ■ ■ {Uuj2)a2ai{Uuji)aiao (A 10) 



(Eao - Ea„_, + Er=l' ^^)iEao - i?a„_. + E7=l ^0 ■ • • (Eao - Ea, + C^i) ' 

In order to describe the line profile for the multiphoton emission, we have to consider first the amplitude for the multiphoton 
emission in the resonance approximation, which can be defined as 

jjem {^ui)aoai(U^^')aia2 ■ ■ ■ {U^^)a,^_ia (All) 

~ {Ea„ - Ea, + UJl)iEa„ - + C^l + W2) • ■ • [Eaa - E^^ + ^^^^ W,) ' 

The resonance approximation for this n-photon emission process assumes actually the existence of the cascade transition a ^■ 



a„_i — )• a„_2 fli ^ flQ. The problem of cascades will be investigated below. An expression, similar to Eq.( A. 1 1 

be defined for the corresponding absorption amplitude 



{Ea, - Ea„ + Ztl ^^)iEao " Ea„_, + Y.Zl ■ ■ ■ (Eao ~ Ea, + C^l^ 



Tj-ab _ y^i^n/aan-i ■ ■ ■ \^u:2ja2ai\^uiijaiao ("AI?"! 



The factors in denominators of Eq.( A.l l| i generate simple poles at the resonance frequencies. These singularities are removed 



(from the real axis) by inserting radiative corrections into the central and all upper electron propagators, i.e. by introducing 
"radiative dressed" propagators. The insertion of the lowest-order radiative corrections (the vacuum polarization, represented by 
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the Uehling potential, and the electron self-energy) in the central propagator (see Fig. 2) yields 



E 



Sl,S2,...S„, 



(Ea, - Es, + OJl){Ea„ - Es, +L01+LU2)--- (Ea, - + J^tl ^0 



(A.13) 



where is the Uehling potential and represents the matrix element of the energy-dependent self-energy operator 



27r 



(-^(|^|))aiKna2 

^-n-En{i-iO) ■ 



Here the shorthand notation 



is used with 



ri2 



(A. 14) 



(A. 15) 



(A. 16) 



ri2 = \ri — \ and the metric tensor .g^^^a- 
In the resonance case Eq.( AT3]i reduces to 



U 



sc,(2n+2) 



(^cJi )aoai )aia-2 ■ ■ ■ (f^Lj„ )an-lti 



(Ko - Ea, + Wl)(£^ao - Ea^ + CJi + ^3) . . . (i^ao - ^^a + ELl 



aan-i • ■ • {Uuj2^a2ai iJ^uJi )aiao 
{Ea, - i?a„_2 + ^0 ■ • ■ {E^o - ^^a, + ^l) 



(A.17) 



Ressumation of an infinite sequence of radiative insertions to all orders of the perturbation theory (geometric progression) 
leads to the following expression for the emission amplitude 



jj-eTn 



{Ea, - Ea, + UJl){Ea, - Ea, + Wi + W2) • - - (-Bao " K + Etl " ^(ELl ^i)) 

with the (in general complex- valued) energy corrections 



U 
aa 



(A.18) 



(A.19) 



i=l 



Expanding the expression for the matrix element of the operator E i nto Ta ylor series around the resonance energy Eag 
Er=i ^» ^ Ea and retaining only the leading term in the correction Eq.( A.19 1 yield 



aoa.1 (f^ai2)°i°2 ■ • ■ (^ti;„ )<in-ia 



where 



(Ko - Ea, + LOl){Ea„ - Ea, + + iJ^) ■ ■ ■ (Ea,, + ELl ^^ ~ K) 
Va=Ea + i^iEa))aa + V^a ■ 



(A.20) 



(A.21) 
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Now let us turn to the insertion in the first upper electron propagator. After performing time and frequency integrations the 
corresponding 5-matrix element reads 



U 



sc. (271+2) _ 



E 



si,S2,...s„, i^ao ~ -S'si + Wi)(i?ao " -^S2 + '^l + ^2) ■ ■ ■ {Eag ~ -£'s„-i + J2i=l 

S7i + l,---,S2n-l,S2n 



(^tJ„ )snS„ + l (t^W„+i )s„ + iS„^.2 • ■ • (^"2n-l )s27i-lS2n ( ^1^2,1 ) 



S2na0 



The resonant case is characterized by the conditions s„_i = s„ = a„_i,si = ai,S2 = a2,...,s„+i = a, s„+2 
a„_i, . . . , S2n = ai, which imply the corresponding scattering amplitude 

(^LJi )cioai (^W2 )'^i''2 • • ■ 1,i-2a„_l 
- + ^l){Ea„ - Ea, + + • • • (Sa„ - i^a„_i + E^l' <^^) 

an-io-ij-^^n )aaji_i • • ■ (^cl;2 )a2ai i^^i )aiao 



(A.22) 



jjres.(2n+2) _ 



(A.23) 



We can assume, that all necessary resonant insertions into the central electron propagator in Fig. 2 have been akeady introduced. 
Repeating the radiative insertions in the upper electron line in the resonance approximation and performing resummation of the 
resulting geometrical progression finally we find 



together with 



n—l n—1 



(A.24) 



(A.25) 



Expanding again the operator S in Eq.( A.25 i into a Taylor series close to the point of the resonance, replacing then Eq.( A.25 1 
by 



and using also Eq.([28]t, we obtain 

jjem 



(A.26) 



(A.27) 



Inserting radiative corrections into the remaining electron propagators and repeating the procedure described above, we finally 
arrive at the following expression for the emission amplitude 



jjem 



aoai {U^^^)aia2 ■ ■ ■ (^tj„)a„_ia 



{Ea„ + C^i - KJ . . . {Ea, + ElTl' - ya^^^){Ea, - E^ + ELl ^^ ^ ^a) 

where energies Va, Va^,..., Va^ are improved by the radiative corrections: 



(A.28) 



Va — Ea + La ~ -Fq , 
Vai -^ai ^" Lai 2"^*^1 ' 



Va^-Ea,^+La„--Ta^, 



(A.29) 
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La, La-^,..., La^ are the Lamb shifts of the levels a, ai,..., an and Fq, Ta^,---, are the corresponding widths. 

As the next step towards the evaluation of transition probabilities one has to take the square modulus of t/"^™, to inte- 
grate over all emission directions v and to sum over the polarization e of each photon. Defining the partial width Ta^a^ of the 
level ai associated with the transition oi — > 02 as 

Ta.a. = / d'^m.aj? , (A.30) 

e 

where w''^'* is the resonant frequency, integrating over all photon directions and summing over polarizations we obtain an 
expression for the 71-photon cascade transition probability 



/ 1 \ 3 r r 



\2^J \Ea, - Ea, + UJi - La, + ^TaA^ \Ea„ - Ea, + UJl + OJ2 - La, + ^TaJ^ 

r r 

a„_2a„-i a„_ia ( A 3 1) 

\Ea„ - Ea^_, + YT'l ^^ " ia„_, + |r,„_, P \Ea, ~ Ea + YJU " La + \T aV ' 



where Fagai,---, ^ and Ta„_ia are the partial widths, associated with the transition ai — )■ ao, a2 = ai,. . ., a„_i ^• 

a„_2 and a — > a„_i. This is the line profile for the cascades transition, when all photons are assumed to be registered. 
From the Eq.( A.31 1 we can obtain also an expression for single-photon transition ai — > in the case when ai,. . . , a„_i, a 
states are excited (unstable) states. For this purpose we perform first integration over uji with fixed 0^2, W3,..., and after this 
perform the integration over uj2 with fixed ^3, a;4,..., a;„ etc. These integrations can be extended along the entire real axis, since 
only the pole terms are contributing. We close the contour in the upper half-plane. The poles for the integration over wi are 
located at 



l-^l ^ — Eai — Eaa + ^ai + 7:^ a. 



'^P^ = Ea2 - Eao - ^2 + La, + -T a. 



^"^ - Ea-Ea,-J2uj^+La+ '-T a (A.32) 



1=2 



2 



Employing Cauchy's theorem and collecting all residue contributions we repeating this procedure for all remaining frequencies 
until LUn- After some algebraic transformation we finally arrive at the expression 



dWa 



2tt 



F F F 



{LJ-^aa„-,r+i(Ta+ra„_,)^ 



(A.33) 



In Eq.(A.33i any explicit dependence on the states ai, ..,a„_2,ao has been disappeared. These states enter only implicitly 
through the definition of the one-photon decay widths. If all states a, ai,..., a„_2 decay via single channels i.e. Taai — Fa, 



^aia; — ctc, the formula Eq.( A.33 1 reduces back to Eq.(3) in the text 



Appendix B. Reemission of the photon, emitted by one atom, by another atom: QED derivation 



We start with the derivation of the Lorentz profile for the emission line from the QED description of the photon scattering on 
the atomic electron 1 13J . lil9J . 

The Feynman graph for the one-photon scattering process is depicted in Fig. 1 for n = 1. The resonance approximation is 



defined by the condition following from Eq. (A. 7 1: 



UJi — UJ2 — Ea — Eao ■ 

A scattering amplitude C/"^ *^^-' corresponding to Eq. ( |A.8| with n = 1 looks like 

siao 



jjsc.(2) ^ 



Eao - Es, 



(B.l) 



(B.2) 
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In the resonance approximation 



and the emission amplitude is defined as 



jjsc.{2) _ {U)aoai{U )aiao 
Eao ^ + W 

jjem.(2) ^ {U*)aiao ^ ^-g 

Thus the emission amplitude can be obtained from the scattering amplitude in the resonance approximation by omitting absorp- 
tion matrix element {U)aasi- 

A transition probability for one-photon emission process a — > ag + 7 follows from the Eq ( |A.33| l and looks like 

dWaao^^- ^""^t ■ (B.5) 

It is assumed that there are no other decay channels for the transition a — > oq apart from the one-photon decay a — >^ qq + 7 
Actually the unique example of such situation is the Lyman-alpha transition 2p — > Is + 7. The normalization condition Eq. (5) 
is valid for the Lorentz profile Eq. (B.5 1. 

Now we will describe the situation when the Lyman-alpha photon, absorbed and emitted by one atom, is then reabsorbed and 
reemitted by another atom. The corresponding Feynman graph is depicted in Fig. 3. The resonance condition Eq. (Bl) is now 
modified as 

J =UO = Ea' - Ea-^ ^Ea~ Ea, , (B.6) 

where quantities with or without "prime" index correspond to the electrons in two different atoms. The S -matrix element, 
corresponding to the graph Fig. 3 looks like 



In this matrix element the variables xi, X2 correspond to one atom and the variables 2:3, X4 correspond to another one. D*^^ 
denotes the transverse photon propagator in the Coulomb gauge. We employ this propagator since we want to describe the 
emission of the real (transverse) photon by one atom and the absorption of this photon by another atom. This propagator can be 
presented in the form Ii4l1 

1 f°° 

D%ixi,X2)^— dr!/;,(|l]|,ri2)e-'"(*i-*^) , (B.8) 



Performing the standard integration over the time and frequency variables and using the relation ( A. 5 1 we arrive as the following 
expression for the scattering amplitude: 

^, {E„, - Ea' + uj'){En - Ea„ + u:) 

In the resonance approximation 

' {Ea- -Ea'„+u')iEa-Eao+Uj) ' ^ ' ' 



According to Eqs. ( |B.3|l-(|B.5| l we obtain the emission amplitude in the resonance approximation by omitting the absorption 

ID 



matrix element in Eq. (B. 

{^e^^-^-^^^y^^ _ (V.a.)(V3a3)^ -'7;:!;°;p-^ ) , , {a2A;i''%ao 

(B.12) 



{Ea'-Ea'+ij'){Ea-Eao+u:) 
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Singularities in the denominators in Eq. (BTTTJi should be avoided by the summation of the electron self-energy and vacuum 
polarization insertions in all the electron lines in Fig. 3. Then, after taking the square modulus of Eq. ( |B.12[ ), summation over 
the polarizations, integration over the angles for emitted photon and finally, after the integration over the frequency lo' we obtain 
the following result (setting Oq = oo = Is, a' = a = 2p) 



27r 



(a; - a;2pis)^ + 



1^1 



s'2p'2pls| , 



(B.13) 



where 



h 



s'2p'2pls 



V ^23 



xVi.(r5)?/'2p(r5) . (B.14) 



In Eq. ( B.14[i the one-electron Dirac wave functions ip2p, V'ls for the electrons in the two different atoms are present. Thus, the 
integral (B.14i depends on the distance between two atoms, i.e. on the density of the atomic gas. It is convenient to fix the origin 
of the coordinate system at the nucleus of an atom which absorbs the photon. 
Then we can present in the form 



f23 



F2 - rsl 



\R- 



(B.15) 



where r', r are the distances between the electrons and the nuclei in the emitting and absorbing atoms, respectively and R is the 
distance between the two nuclei. Assuming that r',r ^ R we replace the distance r23 by R. 
We employ equalities 



{{V 2a2) {"^ Sas) f {r 23)) A' B'AB -{[h2[h3f {r 23)]]) a'^a 
= {EagEa' - EagEa - E^'^Ea' + E ^'^E a) { f {r 12)) a'^a 



Qa'aao ' 



where h{r) is the one-electron Dirac Hamiltonian for the bound electron with arbitrary potential V{r) 

h(r) — — zaV + j3m — eV{r) , 



(B.16) 



(B.17) 



d, /3 are the Dirac matrices, m, e are the mass and the charge of the electron, respectively. The wave functions in the matrix 
elements are assumed to be the eigenfunctions of the Hamiltonian Eq. (B.17i with the eigenvalues Eag, Ea'^, Ea, E'^. In case 
when Ea'g = Eag, Ea' — Ea Eq. (B.16 1 reduces to 



{{^2a2){^3a3)f{r23))a'ga'aag {Eag Ea^ f ij' 23^ a'ga' aag 



Now, employing Eq. (B.18 1 with the evident approximation 

1 

''23 



1 



(B.18) 



(B.19) 



J{E^-E^„)r23 ^pt{E^-E^^)R 



(B.20) 



and taking into account the orthogonality of the wave functions ^pag, i'a, we find that the second term in the brackets in Eq. 
( |B.14| i tm^ns to zero. Then 



Ils'2p'2pls{R) — 77l(<52pls)P 

Jx 



(B.21) 



In Eq. (B.21 1 one of the matrix elements {d)2pis originates from the emitting atom and another matrix element originates from 
the absorbing atom. In the nonrelativistic limit BTl 



\{d2pls)? ^ {E2p - E^,f\{r)2pls\^ 



(B.22) 



and 



2p 



-^{E2p-Ei,f\{r)2pis[' 



(B.23) 
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Hence 

3 r 

We can average the result over the positions of the emitting atoms, surrounding the absorbing atom, assuming the distribution of 
these atoms spherically symmetrical and introducing the density of the emitting atoms p{R). Then 

where 



r . . _ 97r f p[R)n-'dn T% 



is the dimensionless function which represents the frequency distribution for the reemitted photon in Eq. (8). In a simple model, 
employed in section II B the normalization of the function p{R) was chosen as 



T I mp-EM- - 4{E,p-Euy J - ' • ^^-"'^ 
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Fig. 1. Feynman graph describing tlie elastic scattering of n photons on an atomic electron in the state oo 
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Fig. 2. Feynman graph corresponding to lowest-order radiative insertions into the central electron propagator in Fig. 1. 



n = a 



Fig. 3. The Feynman graph describing the resonance photon rescattering. The two vertical solid lines describe the electrons in two 
one-electron atoms. The ground and excited states of the one atom are denoted as Oq, a', the states of another atom oq, a. The internal wavy 
line denotes the transverse photon propagation between two atoms. The Coulomb gauge for the photon propagators is assumed. 



